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EXPERIMENTAL AND ANALYTICAL STUDIES OF CRYOGENIC 

PROPELLANT TANK PRESSURANT REQUIREMENTS 

By M. E. Nein and J. F. Thompson* 

George C. Marshall Space Flight Center 


SUMMARY 


The extensive requirement for pressurization of cryogenic propel- 
lant tanks of launch and space vehicles has directed attention to the 
need for accurate methods of analysis of propellant tank thermodynamics. 
This paper presents the results of experimental and analytical studies 
of pressurization gas requirements for cryogenic liquids. Experimental 
results are analyzed for cylindrical and spheroidal tanks ranging in 
size over four orders of magnitude. A parameter study of the con- 
trollable variables of a pressurization system design illustrates their 
effect on ullage gas temperature. 


Pressurization data are provided for use in the development or 
checkout of analytical pressurization models and for design of pressuri- 
zation systems for future launch and space vehicles. A tank pressuri- 
zation computer program, utilizing recommended coefficients, can be 
used to predict total and transient pressurant requirements andjjlitfge" 
temperature gradients within 10% accuracy, 

INTRODUCTION 



Determination of the pressurant gas weight for cryogenic propellant 
tanks is complex and defies exact analytical treatment because of the 
interdependent transient phenomena of heat and mass transfer that occur 
simultaneously in a propellant tank. Mathematical models describing 
the internal thermodynamics of tank pressurization have been developed 
by various investigators. 


The analysis by Clark (Ref. 1) represents an analytical solution of 
the governing equations that predict the transient temperature, the 
response of the pressurant gas, and container wall. However, the 
solution requires assumptions, such as constant tank pressure and zero 
initial ullage, that are not always met with real systems. The studies by 
Coxe and Tatum (Ref. 2) are based on analysis of a system in which the 


^Assistant Professor, Mississippi State University, Aeronautical 
Engineering Department. Formerly with Propulsion Division, MSFC. 


ullage gas is thermally mixed and heat transfer between the gas and 
the wall is independent of time and space. Gluck and Kline (Ref. 3) 
used dimensional analysis to express gas requirements as a function 
of known system parameters; they determined, experimentally, 
quantities of interfacial mass transfer and gas phase heat transfer. 

Epstein (Ref. 4) presented a numerical method for calculation of 
pressurant gas requirements that contains a number of phenomena 
absent from previous analytical methods. However, empirical data 
are required to evaluate many constants and physical parameters. 

To provide a reliable method for determination of pressurant gas 
requirements, the experimental data on pressurization obtained by the 
Marshall Space Flight Center during the SATURN launch vehicle 
development were applied to the method of Epstein. The physical para- 
meters and the previously indeterminate constants were developed. 
After modification, this numerical method is capable of accurate pre- 
diction of pressurization gas requirements for cryogenic propellant 
tanks. 
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PRESSURIZATION REQUIREMENTS AND LAUNCH VEHICLE DESIGN 


The increasing size and complexity of space launch vehicles 
necessitates optimization studies of the many subsystems involved in 
launch vehicle design. The propellant tank pressurization system is 
of particular importance because its weight is large in comparison to 
the weight of other subsystems. Weight optimization studies of pro- 
pellant tank pressurization systems for the SATURN V S-IC stage were 
used to establish the location of the oxidizer and fuel tanks within the 
overall vehicle configuration {FIG 1 and Ref. 5). Even the pressuriza- 
tion system components such as heat exchangers, pressurant lines and 
controls, weigh considerably less than the pressurizing gas. 

A further indication of the need for optimization of pressurant 
requirements is illustrated in FIG 2. The pre s surant- ma s s / tank- 
pressure ratios of typical launch vehicles is given as a function of 
vehicle thrust, thrust being representative of vehicle size. Although 
there is a great deal of difference between the propellant tank configura- 
tions of tactical missiles and space launch vehicles, a near linear 
increase occurs in pres surant-ma ss / tank-pr e s sure as vehicle size 
increases. Considering only pressurant gas weight, it appears advan- 
tageous to use helium as a pressurant. If, however, the weight of the 
pressurant storage containers is included in the weight of the pressuri- 
zation system, the use of helium as a pressurant in most instances 
results in a weight penalty. 

For vehicles with high acceleration and low turbo-pump NPSH 
requirements, it is possible to eliminate the pressurization system, 
relying only on the self pressurization of the saturated propellant 
(flash boiling). However, flash boiling pressurization results in high 
pressurant weight and can only be justified if it significantly simplifies 
vehicle design. Because of the infant knowledge of cryogenic tank 
pressurization at the initiation of the Saturn launch vehicle develop- 
ment program, a long series of pressurization experiments was con- 
ducted at MSFC to obtain system design information and scaling laws 
for the large propellant tanks of the Saturn I vehicle. Results of this 
experimental program and correlations with analytical studies are 
presented in the following sections of this report. 
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EXPERIMENTAL PROGRAM 


Test Facilities 

The experimental work was conducted on five tank configurations 
at the Marshall Space Flight Center: 

A. Saturn I, S-I Stage, Multiple Interconnected LOX Tanks 
(FIG 3a) 

B. Saturn I, S-IY Stage (FIG 3b) , LOX and LH 2 Tanks 

C. A 6. 5 x 39 foot (DxL) cylindrical LOX tank (FIG 4) 

D. A 13 x 26 foot (DxL) cylindrical LOX tank (FIG 5a, FIG 5b) 

E. A 1 x 3 foot (DxL) cylindrical LOX tank 

The test parameters for these tank configurations are compared 
in Table I. Configurations A and B were flight vehicles and thus con- 
tained the standard test instrumentation of the Saturn propellant feed 
system, including continuous liquid level sensors, tank pressure, 
pressurant flowrates, and supply temperature measurements. Con- 
figurations C, D, and E were equipped with many thermocouples along 
the tank axis. Thermocouples, mounted at several radii at three 
elevations in these tanks, allowed measurement of radial temperature 
gradients. Wall temperatures were measured in Configurations C 
and D by thermocouples on the inside and outside surfaces of the tank 
at several locations. The locations of the temperature sensors in 
these tanks are shown in FIG 6a, 6b. Special calorimeter plates were 
mounted in both tanks for determination of gas-to-wall heat transfer 
coefficients. Finally, gas sampling devices were placed at several 
locations to measure ullage gas concentration gradients. 

Configurations C, D, and E were equipped with heat exchangers 
that provided a variable pressurant inlet temperature up to 1000°R. 

The pressurant gas was introduced at the top of the container through 
a distributor (either a deflector plate - Configuration C and E, or a screen 
arrangement- Configuration D) to minimize inlet velocities and distur- 
bances of the liquid surface by impinging gas jets. FIG 7 shows a typical 
distributor configuration. Pressurant velocities at the distributor 
periphery are given in Table I for the five test configurations. 
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The tank Configurations C, D and E could be sloshed at rotational 
or translatory oscillation in excess of the first critical frequency of 
the tank. Configurations A, C, and D were equipped with cameras so 
that the conditions inside the tank could be observed. The results of 
tests conducted with the five tank configurations are presented in 
FIG 13-38. The conditions of these tests are summarized in Table II. 

Instrumentation 


Analysis of ullage gas temperature history required a temperature 
probe with fast response characteristics and good accuracy. A fast 
response temperature probe (FIG 8A) was designed at MSFC consisting 
of a fork-like support with a 30 gage CuCo welded thermojunction. The 
length-to-diameter ratio of the thermocouple wire and its distance from 
the fork base were determined using an analog computer representation 
of the heat transfer conditions around the probe assembly. FIG 9A shows 
the response time; 63. 2% of the total temperature change was attained in 
eight seconds when the probe was extracted from liquid oxygen into a gas 
circulating at a velocity of about 3 ft/ sec. The response of the probes 
during a pressurization test was also determined (FIG 9B); the fork-type 
thermocouple has a good response characteristic. 

A thermocouple mounted on a long, rod-like support (FIG 8B), which 
was designed for liquid measurement in the high vibration environment 
of static and flight testing, exhibited an extremely poor response in the 
gas phase as indicated in FIG 9A. Response time to 63. 2% of total temper- 
ature change was in excess of 10 minutes. Commercial temperature 
probes of the resistance thermometer type (FIG 8C) were also investigated 
under these conditions. Although their response was considerably better 
than the flight type thermocouple (63. 2% temperature change in approxi- 
mately 50 sec), it was too slow for the pressurization studies. 

Pressure measurements in the ullage space, pressurant supply lines, 
and liquid discharge lines were made with close-coupled pressure trans- 
ducers to assure good response characteristics. The pressurant flowrate 
and liquid discharge flowrate measurements were obtained with turbine 
type flowmeters. Liquid level before and during the tests was measured 
by capacitance discrete level probes and continuous delta P measurement, 
of the liquid column. 
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Test Results 


a. Heat Transfer Coefficients 

Heat transfer between pressurant and tank side walls was 
measured during pressurization tests in Configuration C by two plate 
calorimeters. Each calorimeter was a 12 x 12 inch, 30 gage copper 
plate mounted from teflon spacers parallel to and at a distance of 
four inches from the tank wall (FIG 10). Three thermocouples, spaced 
to represent equal calorimeter areas and connected as a thermo-pile, 
provided a temperature /time history of the copper plate before and 
during the tests. The local ullage gas temperature was measured in 
the vicinity of the calorimeter (FIG 11). The calorimeters were 
located 1 1 ft and 30 ft from the top of the test tank. 

For determination of heat transfer coefficients, it was assumed 
that heat transfer to the back side of the plate (towards tank wall) was 
by free convection because of the shielding effect of the plate-to-wall 
arrangement. The free convection coefficients for a one component 
gas were evaluated by the equation of Jackson and Eckert (Ref. 7); the 
results are plotted in FIG 12. The free convection heat transfer 
coefficient was also calculated for two component mixtures based on the 
time and space dependent helium-oxygen concentration in the tank. 

The tota.1 heat transfer to the calorimeters was then corrected using the 
calculated free convection effect on the back side. The heat transfer 
coefficients to the front of the calorimeter plates measured in Tests 
130-9, -10, -15 are present in FIG 13 and FIG 14 using gaseous oxygen 
and helium as pressurants. Ullage gas-to-wall heat transfer was also 
evaluated from wall temperature measurements at a location 3. 5 ft 
from the top of the tank. Wall measurements at locations initially 
below the liquid surface produced erroneous readings and were dis- 
carded. These coefficients were corrected by subtracting the effect of 
external heat flux from the measured wall temperature rise. During 
a flash-boiling test, which did not require pressurant flow, the wall 
temperature rise indicated an external heat flux of 1 3 RTU/min ft 2 ; 
this compares very favorably with a calculated flux of 1 5 BTU/min ft 2 
(Ref. 8) and confirms the method used for correcting wall measurements. 


Inspection of FIG 13 and FIG 14 shows very good agreement 
between measured and calculated heat transfer coefficients. It is 
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noted that the gas-to-wall heat transfer coefficient is definitely within 
the forced convection regime for the oxygen tests, but in the free 
convection regime for the helium test. Although the heat transfer 
coefficient by forced convection diminishes with increasing distance 
from the pressurant distributor, the free convection contribution (see 
Equation 1) compensates for this decay to such a degree that a nearly 
constant heat transfer coefficient is obtained along the tank bulkhead 
and side wall. 

b. Sloshing Effects 

Pressurization studies conducted at MSFC have shown that 
there is little benefit derived from the use of helium as a main pressur- 
ant for cryogenic propellants. However, it was determined experi- 
mentally that prepressurization with helium reduces pressure decay 
during liquid sloshing near the critical frequency. It is assumed that 
the helium acts as a buffer zone between the splashing cryogenic liquid 
and the condensable pressurant, suppressing excessive mass transfer. 

FIG 15 shows a typical tank pressure history for a stationary 
liquid oxygen test tank as compared to a pressure history in which the 
liquid sloshes near the first critical mode of oscillation (Ref. 9). The 
tank was prepressurized, with either helium or nitrogen, followed by 
main pressurization during liquid expulsion with super-heated oxygen. 
The tank pressure history during the slosh test (using helium as a 
prepressurant) is nearly identical to the pressure history of the non- 
sloshing expulsion test. In contrast, prepressurization with gaseous 
nitrogen resulted in a marked pressure decay during the sloshing of the 
liquid, which was not evident during a non- sloshing expulsion test with 
gaseous nitrogen prepressurization. 

c. Ullage Gas Concentration Gradients 


Gas flow conditions and the concentration of helium gas in a 
cryogenic propellant tank during pressurization discharge were studied 
in test Configurations C and D. Spectrographic analyses were made of 
gas samples taken at various positions in the tanks. Samples taken at 
various elevations in tank Configuration C just before the end of the 
tests yielded the results shown in FIG 16. In the test in which helium 
was used for prepressurization and oxygen as the main pressurant, 
the helium concentration is maximum at 1 2 ft above the liquid, and 
gradually decreases in both directions. 
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The concentration of oxygen near the liquid surface is probably 
caused by accumulation of the gaseous oxygen that is initially in the 
ullage before prepressurization. For comparison, FIG 16 also shows 
the concentration of helium above the liquid oxygen for the case in which 
helium prepressurization is followed by pressurization with helium 
during liquid expulsion. The oxygen concentration at 10 ft above the 
liquid interface was only 6 percent by volume. The total amount of 
gaseous oxygen in the ullage was only slightly larger than the amount 
of oxygen in the ullage before prepressurization (0. 77 moles vs 0. 73 
moles). This indicates that interfacial mass transfer, although small 
under these conditions, was in the form of evaporation. 

d. Mass Transfer 


A comparison of mass transfer results obtained in Configur- 
ation C with results obtained by Clark (Ref. 1) is shown in FIG 17. 
Condensation in excess of 30 percent of the pressurant flow was found 
by Clark during liquid nitrogen expulsion tests with a 1 x 3 ft cylindrical 
tank. Similar results were obtained with the MSFC test Configuration 
E, also shown in FIG 17. The mass transfer measured in test Con- 
figuration C indicates that condensation was 5-10 percent. Condensation 
in the larger facility is less because of the smaller wall-area/volume 
ratio of a larger tank. 

Comparing the condensation in the small tank with that in the 
large tank on the basis of wall-area/volume ratio, the values are 
approximately equal. During tests at high pressurant inlet temperature, 
initial evaporation noted in Configuration C diminished as the test 
proceeded. However, Clark had found increased condensation at higher 
pressurant inlet temperatures in small tanks. These conflicting results 
point out the incomplete knowledge of mass transfer. 

e. Condition of Liquid Interface 


The condition of the liquid interface in Configuration C and 
during the launch and flight of SA-5, are shown in several frames from 
a movie taken inside these tanks (FIG 18 and 19). Violent boiling 
occurred during venting of the tank before prepressurization. As the 
vents were closed and prepressurization proceeded, the liquid surface 
became nearly quiescent before discharge. After discharge began, 
disturbance of the liquid surface caused by pressurant flow and 
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acceleration of the liquid surface were observed; the disturbance 
diminished as time and distance between the surface and the pressurant 
inlet increased. 

f. Radial Ullage Temperature Gradients 

Radial temperature gradients obtained with Configurations C 
and D are shown in FIG 20. In both cases the radial gradients were 
small, and there apparently exists little difference between the gas 
flow conditions in. the two tanks, even though the gas distributors, 
baffling, and tank diameters are not comparable. 

The temperature probes at X/D-0. 025 in Configuration D, which are 
located between the antislosh baffles (FIG 5b), recorded virtually the 
same temperature as probes at smaller radii. It was concluded that 
the gas circulation in the tank is not appreciably affected by the anti- 
slosh baffles, and subdivision of the tank into volume elements perpen- 
dicular to the tank axis is permis sable for the pressurization analysis. 

g. Axial Ullage Temperature Gradients 


The axial ullage temperature gradients obtained in tests 
130-6 and 130-7 with Configuration C (FIG 21a, b, 22a, b) became 
approximately linear as the test proceeded. These two tests were con- 
ducted with oxygen as pressurant at about 550°R. There was a rapid 
increase in temperature of about 30°R immediately above the liquid 
interface in these tests, indicating that mass transfer was small. In 
tests 130-9 and 130-10 (FIG 23a, b and 24a, b) with the same Configuration 
with oxygen pressurant at a lower temperature, the ullage temperature 
gradients are much flatter; the rapid increase in temperature immedi- 
ately above the liquid interface is still in evidence. The ullage 
temperature gradients in this same configuration witn helium as pres- 
surant (test 130-15; FIG 25a, b) are concave, rather than linear as in 
the tests with oxygen as pressurant, and the increase in temperature 
just above the liquid interface is very gradual. The concave shape is 
to be expected in this case because the mass transfer is in the form of 
evaporation with an ullage that is predominatly helium. The linear 
ullage temperature gradients in the tests with oxygen as pressurant 
indicate that the mass transfer is very small with an ullage that is 
predominately oxygen. 
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h. Other Test Results 


Tests are being performed with Configuration D, but so far 
only three tests have been completed. The pressurant distributor in 
this configuration was designed to minimize the gas circulation in 
the tank, reducing forced convection heat transfer. While this is the 
desired condition for optimum pressurization system operation, it is 
detrimental to the response time of the temperature probes as the 
liquid interface passes. Precise ullage temperature gradients will 
not be available until this instrumentation is improved. However, 
preliminary data, with very hot GOX used as pressurant, indicate that 
the temperature gradients are concave rather than linear as was the 
case in the tests with Configuration C using colder GOX as pressurant. 
The concave temperature gradients found in the helium pressurant 
tests with Configuration C were also in evidence with Configuration D. 
Pressurant flowrates and wall temperature gradients from these tests 
are presented in FIG 26a, b, c, 27a, b, c and 28a, b, c . 

Pressurant flowrates in the LOX tanks of the SATURN I, S-I 
stage, during static test and flight are presented in FIG 29a and 30. 
FIG 31a and 32a show pressurant flowrates in the LOX and LH 2 tanks 
of the SATURN I, S-IV stage, during static test. Finally, ullage 
temperature histories obtained in a very small tank, Configuration E, 
containing LN 2 pressurized with nitrogen are given in FIG 33. 

PRESSURIZATION ANALYSES 


Previous Work 

Pressurized discharge from cryogenic liquid containers was studied 
analytically and experimentally by Clark (Ref. 1) under sponsorship of 
the Army Ballistic Missile Agency and later MSFC. The analytical 
solutions obtained by Clark were applied to test data obtained for Con- 
figuration C. In FIG 34 the axial temperature gradient through the ullage 
gas is shown as a function of distance from the tank top or gas distributor. 
Excellent agreement with test results was obtained for an assumed gas- 
to-wall forced convection heat transfer coefficient of 10 BTU/hr ft 2 °R. 
Agreement for a coefficient of 2 BTU/hr ft 2 °R, approximately in the 
range of free convection, was poor. This illustrates one limitation of 
analytical solutions in which the gas-to-wall heat transfer coefficient 
enters as an independent variable. 
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In spite of this restriction and the assumption of initial zero ullage 
volume, the method by Clark was successfully applied in design analy- 
ses of the Saturn I pressurization system. While Clark’s analysis 
assumed stratification of the ullage gas and constant heat transfer 
coefficient, the analysis by Coxe and Tatom (Ref. 2) was based on the 
assumption of a complete thermally mixed ullage gas and constant heat 
transfer coefficient. FIG 35 and 36 compare test results obtained with 
MSFC Configuration C with the analytical predictions by the method of 
Coxe and Tatom. Toward the end of the test, agreement is good 
possibly because the conditions of constant heat transfer coefficients 
are approached in the large ullage near the end of the run. 

A comparison of the pressurant flow requirements with predictions 
by an analog computer simulation developed by MSFC, is shown in 
FIG 37. Representation of the pressurization thermodynamics by ana- 
log method was difficult because of scaling problems and the extreme 
sensitivity of the equations to tank pressure fluctuation. In FIG 35, 36, 
and 37 pressurant flow requirements are also compared with a digital 
computer program developed by Rocketdyne (Ref. 4) and modified by 
MSFC (Ref. 10). This program closely matches test data. However, 
the program insufficiently describes mass transfer and is sensitive to 
fluctuations of ullage pressure. These fluctuations do not appear in the 
measured flowrates because they are apparently counteracted by the 
effects of evaporation and condensation (Ref. 11). 

Summary of Rocketdyne Program 


This program makes maximum use of the techniques of digital com- 
puter calculations and is not subject to the restrictive assumptions that 
are made in other programs. This method was therefore chosen by 
MSFC for pressurization system analyses. However, extensive com- 
parisons of the program with test data were required to evaluate the 
physical parameters and constants initially contained in the program as 
indeterminate identities. The equations were modified when necessary. 

This program includes in its calculations a pressurant gas storage 
tank, heat exchanger, and flow control valve. It considers a propellant 
tank with or without outside insulation and pressurized with either 
evaporated propellant or with a gas stored under pressure in a storage 
tank in which the gas expands nonadiabatically. The ullage pressure is 
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controlled by a pressurant flow control valve that has finite maximum 
and minimum areas and may be either the on-off or the continuously 
regulating type. In the propellant tank the ullage gas may be a two 
component mixture of evaporated propellant and another gas. The 
ullage gas temperature, composition, and properties are considered 
functions of time and of axial, but not radial or circumferential, 
distance. Liquid and wall temperature and properties are treated in 
the same manner. The heat transfer modes considered are shown in 
FIG 38. Mass transfer within the ullage and at the gas-liquid interface 
is considered. The effects on heat and mass transfer caused by gas 
circulation, as influenced by pressurant gas inlet velocity, is also taken 
into account. 

Modifications in the Program 


In the course of the comparisons with test data, it was necessary 
to make several modifications in the program to obtain good data 
correlations. These modifications are discussed in Reference 10. 

The ullage gas-to-wall heat transfer coefficient, which decreases 
exponentially from the tank top, is written as the sum of a free con- 
vection coefficient and a forced convection coefficient.'* 5 

/ Vd \ ■ PwZ 

h gw = h c + h o \v od y e < X > 

where h Q is an input constant. 

Thus the forced convection coefficient at the tank top is a linear 
function of the pressurant volumetric flowrate (V^)from the distributor. 
The free convection coefficient(h c )is calculated by the free convection 
equation, 


h c = c . k 


g c c 7 CpP (|Tg_-T w M| 


L w' 


T g M- k 


(2) 


^Schmidt (Ref. 12) also writes the total heat transfer coefficient as 
the sum of the free and forced convection coefficients. 
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In the same manner the gas -to-liquid heat transfer coefficient at 
the gas -liquid interface is written _ y 

~ P S 1 

' e (3) 


^s ^sc h so 


4^Y 


where h so is an input constant. 


The free convection coefficient h is calculated by the equation 


h = c k 

SC 1 


gc c 7 c P P 2 (lT g -T L | ) 


T g ^k 


(4) 


It was found that both forced convection coefficients at the tank 
top could be calculated more accurately by a forced convection equation 
of the standard form expressing the Nusselt number as a function of the 
Reynolds and Prandtl numbers: 



Thus, the ullage gas-to-wall heat transfer coefficient and the gas-to- 
liquid heat transfer coefficient at the gas-liquid interface are better 
calculated according to equations 7 and 8. 



~Pw Z 


hg W 

= h c + h D e 

(V) 


- Ps Zi 

(8) 

h s 

- h gc + h so e 


where h Q and h go are calculated by equations 5 and 6, rather than 
being input as constants, and h c and h sc are calculated by equations 2 
and 4 . 
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It was also found that the liquid-to-wall heat transfer coefficients 
could be better calculated according to a free convection equation 
rather being taken as constant: 


c 4 


c , k 


g c 7 (3 c p z ( T - T ) 
6 c 7 k LP p L w 


H k 


(9) 


As in the case of gas-to-liquid heat transfer coefficient at the 
gas-liquid interface, the mass transfer coefficient at the interface 
was written 


Y = Y + Y 
1 s sc so 


(ty 


- Ps^i 


where Y so is an input constant. 


(10) 


The free convection coefficient (Y sc ) is calculated by the equation 


^sc c i^ 


g C c 7 P (f^g - <^i 1 ) 


D 


(ID 


The forced convection mass transfer coefficient at the tank top can 
be better calculated by a forced convection equation expressing the 
Sherwood number as a function of the Reynolds and Schmidt numbers: 


Y S o r 

D 



( 12 ) 


Thus, the mass transfer coefficient at the gas-liquid interface is 
calculated by 


- p s z. 

Ys ~ Y sc + Y so e (13) 

where Y so is calculated according to equation 12, rather than being 
input as a constant, and Y sc is calculated by equation 11. 
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Evaluation of Program Parameters 


All pressurization analyses contain numerous parameters that 
must be known before pressurization requirements can be predicted. 
These parameters determine the heat and mass transfer coefficients 
and the distribution of these coefficients over the tank. Therefore, 
studies were conducted to determine the relative importance of each of 
the parameters involved in the Rocketdyne program, and extensive 
comparisons with the results of the tests were made to obtain numerical 
values for these parameters. 

A summary of the test conditions is given in Table II, and the 
values of the important parameters are given in Table III. The expo- 
nential decay coefficients ]3 w and 3 s in equations 7, 8 and 13 are scaled 
by the equation: 

p = 0.00117 r 2 (14) 

The parameters not listed in this table are of small importance and 
may be taken as zero. 

Comparison with Test Data 

The pressurant flowrate and ullage and wall temperature gradients 
predicted by the computer program using the calculated constants 
from Table III are compared with test data (Ref. 13, 14, 15, 16 and 17) 
in FIG 21-30. In all comparisons the ullage pressure, liquid drain 
rate, ambient heat transfer coefficients, and ambient temperature 
were input to the computer as functions of time. Either the pressurant 
inlet temperature or the heat exchanger performance curve was also 
input. 

FIG 21-25 show comparisons with test data obtained with Configur- 
ation C described in Table I and shown in FIG 4. As can be seen from 
these figures, the agreement between the computer predictions and the 
test data is generally good. The irregularities in the computed pressur- 
ant flowrate, particularly marked in tests 130-6 and 130-7 (FIG 21 and 
22), are caused by the over -sensitivity of the program to changes in 
the slope of the ullage pressure curve. Both ullage pressure curves 
of tests 130-6 and 130-7 have depressions in the latter half of the runs, 
while the slopes of the ullage pressure curves of the other tests were 
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nearly constant. The agreement between the computed and measured 
ullage temperature gradients was good throughout the run for all the 
tests using oxygen as pressurant. In the test with helium as pressurant 
130-15 (FIG 25) the pressurant flowmeter failed. Storage bottle 
pressure and temperature history were used for calculation of an 
average flowrate. Therefore, it was not unexpected that the computed 
flowrate was somewhat below this value. However, the agreement 
between computed and measured ullage temperature gradients was not 
as good in this test as in the test with oxygen as pressurant. This was 
probably caused by deficiencies in the program’s mass transfer calcula- 
tions from the assumption that all heat transfer from the ambient to the 
propellant is converted to sensible heat rather than latent heat. In 
tests 130-9 and 130-10 the ullage heat transfer coefficients were 
calculated from calorimeter measurements and were compared with 
those calculated by the computer. Although the assumption of exponential 
decay of the ullage heat transfer coefficient with distance from the tank 
top (Equation 7) seems arbitrary, the results were in excellent agree- 
ment with the measured heat transfer coefficients (FIG 13 and 14). 

In comparing the velocity decay of a free jet (FIG 39, discussed 
in Reference 18) it was found that the exponential decay of the forced 
convection heat transfer coefficient expressed as a velocity decay 
( v z/v Q ) * is bracketed by the velocity decay of a free jet discharging from 
a circular opening and that of a free jet discharging from an infinite slit. 
This is analogous to the pressurant entering the tank through the gas 
distributor. 

Comparisons with data from the LOX tanks of the Saturn I, S-I 
stage during static test and flight are presented in FIG 29 and 30. The 
agreement between computed and measured pressurant flowrate and 
pressurant inlet temperature is excellent. Ullage temperature measure- 
ments were not available in these tests because instrumentation on flight 
vehicles is limited. FIG 29 shows a comparison of the computed and 
measured flowrate from the flight of SA-5. The agreement was generally 
good, though not as good as in the static test of SA-6. Evaluation of 
SA-5 pressurant requirements was complicated by the complex air flow 
pattern around and between the propellant tanks of the Saturn I, S-I stage 
during flight. The aerodynamic heating was difficult to evaluate; the only 
possible approach was to use average values for all propellant tanks. 

FIG 31 and 32 show comparisons with data from the LOX and LH z tanks 
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of the Saturn I, S-IV stage during static test. These tanks are not of 
ordinary cylindrical shape, as can be seen in Table I; the LOX tank is 
an oblate spheroid and the LH 2 tank contains a convex inward lower 
bulkhead. By computer variation of the characteristic tank radius 
used in equations 5, 6, and 12, it was determined that the proper 
characteristic value should be about two- thirds of the maximum radius 
for the LOX tank. This assumption is theoretically justified because 
a cylinder having the same volume and surface area as an oblate spher- 
oid has a radius equal to 0. 63 times the maximum radius of the oblate 
spheroid. The agreement between computed and measured pressurant 
flowrate in the LOX tank is excellent, as shown in FIG 31. Since the 
pressurant flowrate in the LH 2 tank was a step function, it could not be 
matched at all times. However, the general range of flowrate, as 
computed and measured, is the same, and there is excellent agreement 
between the computed and measured total pressurant weight. 

Test results with Configuration D are shown in FIG 26-28. This 
tank is an approximate 1/3 scale model of the Saturn V, S-IC stage, 

LOX tank. It is the largest single cylindrical LOX tank from which 
test data is currently available. Comparisons of the computer predic- 
tions with data obtained from three tests with this configuration is good 
for pressurant flowrates and tank wall temperatures. 

The final comparison presented is with data from a very small ' 
cylindrical tank (one foot in diameter and three feet long) with flat 
bulkheads (Configuration E). Although pressurant flowrate measure- 
ments were not available in this test, the computed and measured ullage 
temperature histories are compared in FIG 33. The agreement is not 
as good as obtained in Configuration C, probably because equation 14 
for the scaling of the exponential decay coefficients was developed for 
tanks with rounded rather than flat bulkheads. 

Conclusions from Comparisons with Test Data 


These comparisons with test data cover a range of conditions, 
using oxygen, helium,and nitrogen as pressurants and liquid oxygen, 
liquid hydrogen, and liquid nitrogen as propellants in tanks ranging in 
size over four orders of magnitude. The tank shapes were representa- 
tive of those commonly used in space vehicles, namely cylinders with 
various bulkhead shapes and oblate spheroids. As a result of the 
evaluation of the many physical parameters and constants involved in 
the equations, this program can be used to predict total and transient 
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pressurant flow requirements, ullage and wall temperature gradients, 
and gas-to-wall heat transfer coefficients with an accuracy of ± 5%. 

The numerical values of parameters recommended by MSFC for use 
in the program are given in Table III. There are presently no other 
values available in the literature. The characteristic dimension 
used in the calculation of the exponential decay coefficients was taken 
as the radius of the cylindrical section for cylindrical tanks. For tanks 
of other shapes, some comparison with test data was necessary to 
determine the proper choice for the characteristic radius. A value of 
two-thirds of the maximum radius appears acceptable for oblate 
spheroids. 

The comparison with test data indicates a sensitivity of the program 
to sudden changes in ullage pressure. However, in most cases vehicle 
design pressures are either constant or vary in a monotonic manner. 

It was further found that considerable experimental experience with 
pressurization systems is required before this method of analysis can 
be applied reliably to evaluate a new system. 

THE EFFECTS OF SYSTEM PARAMETERS 
ON PRESSURANT REQUIREMENT 

Weight optimization of propellant tank pressurization systems 
demands that a low pressurant density be maintained in the ullage space; 
this is analogous to using a gas of low molecular weight and maintaining 
a high ullage mean temperature. Therefore, 30 pressurization tests 
and 120 computer predictions were used to separate the relative signifi- 
cance of various controllable parameters of pressurization systems and 
to determine their influence on mean ullage temperature. FIG 40 
presents a graphical illustration of the relative influence of these para- 
meters. 

From a central origin, representing a reference condition (Saturn V, 
S-IC Stage) for all parameters, the increase (+ Y) and decrease (- Y), of 
the ullage mean temperature at cutoff is shown as a function of variation 
of the parameters on the abscissa. The parameters were varied over 
a range expected for vehicle design. Thus, pressurant inlet tempera- 
ture can increase or decrease by a factor of two from the reference 
condition, pressure by a factor of three, tank radius by a factor of two, 
expulsion time by a factor three, etc. It was indicated that the pres- 
surant inlet temperature exerts the greatest influence on the ullage 
mean temperature. Diminishing return of this effect did not exist 
within the range of investigation (530°R to 1200°R). The mean temper- 
ature increased as the ullage pressure was increased and also as the 
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tank radius was increased. Increasing the tank wall thickness, heat 
capacity, or density caused a decrease in the mean temperature. 

The pressurant distributor flow area (A D ) that controls the gas-to- 
wall forced convection heat transfer coefficient had a significant 
effect on the mean temperature when the area was reduced, but no effect 
at all when flow area was increased. This indicates that the pres- 
surant inlet velocity for the reference systems was chosen at an 
optimum point. Figure 40 also indicates that helium pressurant must 
be introduced into a tank at a temperature 1. 1 times higher than 
oxygen pressurant to obtain the same ullage mean temperature. This 
confirms the results of other studies (FIG 2) indicating that the 
benefits derived from a helium pressurization system are not based 
on weight optimization. 

CONCLUSIONS AND RECOMMENDATIONS 

a. Pressurization data from cylindrical and spheroidal tanks 
ranging in size over four orders of magnitude are available for 
development or checkout of analytical pressurization models and for 
design of pressurization systems for future launch and space vehicles. 

b. The Rocketdyne tank pressurization program, modified as 
described herein and utilizing recommended coefficients, can be 
used to predict total and transient pressurant requirements and 
ullage temperature gradients with an accuracy of ± 5%. 

c. No significant radial ullage temperature gradient occurs, even 
in tanks with anti- slosh baffles. This permits the assumption of one- 
dimensional stratification of the ullage gas for analytical representation 
of pressurant requirements. 

d. Heat transfer between pressurant and tank walls can differ 
significantly from free convection, depending on tank geometry and 
distributor design. 

e. The strongest influence on pressurant weight is exerted by 
pressurant inlet temperature, for which no diminishing return occurs 
within a temperature range compatible with tank materials. Other 
important influencing factors are tank radius, distributor flow area, 
expulsion time and aerodynamic heating. The effect of wall heat 
capacity is not as significant as might be expected. 
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f. Mass transfer for large tanks is less than previously reported. 

g. Additional work is necessary to develop better techniques for 
measuring gas concentration gradients and mass transfer. 
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TANK CONFIGURATIONS AND TEST PARAMETERS 

































TEST FACILITY ULLAGE PRESSURE (pslo) INLET TEMPERATURE PR) PRESSURANT PRE- 
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2) Sloshing during SA-5 flight unknown 










































































FIGURE I. COMPARISON OF THE WEIGHTS OF THE PROPELLANT 
FEED SYSTEMS OF TWO FLIGHT VEHICLES 
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FIGURE 3B 


SATURN I, SHV STAGE 


FIGURE 3A 


SATURN I, S-l STAGE 
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FIGURE 5a. INTERIOR OF TANK CONFIGURATION D 
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FIGURE 5b. TEST FACILITY FOR TANK CONFIGURATION D 
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FIGURE 9. COMPARISON OF TEMPERATURE SENSOR RESPONSE TIME 
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FIGURE 11. TEMPERATURE RESPONSE OF COPPER PLATE CALORIMETER 
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AND GN 2 PRE-PRESSURANTS UNDER LIQUID SLOSH AND NON 
SLOSH CONDITIONS IN TANK CONFIGURATION C 
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FIGURE 16. MEASURED AND COMPUTED ULLAGE GAS CONCENTRATION 
GRADIENTS IN TANK CONFIGURATION C 




(LB/LB) VS TIME T (SEC) 







■u u M 


mm ■ ■ 


CO. J 




d> IF? 


jr 

ZM 

> i 

i mtr Jm 


45 


FIGURE 18. LIQUID SURFACE CONDITIONS DURING PRESSURIZATION TEST 
IN TANK CONFIGURATION C 
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FIGURE 19. LIQUID SURFACE 
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FIGURE 20. EXPERIMENTALLY DETERMINED RADIAL TEMPERATURE GRADIENTS 
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FIGURE 21c. COMPARISON BETWEEN EXPERIMENTAL AND COMPUTED PRESSURANT 
FLOWRATE, TANK CONFIGURATION C, TEST 130-6, OXYGEN AS 
PRESSURANT 
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TEST 130-9, OXYGEN AS PRESSURANT 
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FIGURE 24d. ULLAGE PRESSURE AND PRESSURANT INLET HISTORIES, TANK 
CONFIGURATION C, TEST 130-10, OXYGEN AS PRESSURANT 
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FIGURE 25d. ULLAGE PRESSURE AND PRESSURANT INLET TEMPERATURE 
HISTORIES, TANK CONFIGURATION C, TEST 130-15, HELIUM 
AS PRESSURANT 





TEST 



FIGURE 26a. COMPARISON BETWEEN EXPERIMENTAL AND COMPUTED 
PRESSURANT FLOWRATE, TANK CONFIGURATION D, TEST 
C 003-7a, OXYGEN AS PRESSURANT 
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FIGURE 27a. COMPARISON BETWEEN EXPERIMENTAL AND COMPUTED 
PRESSURANT FLOWRATE, TANK CONFIGURATION D, 

TEST COO 3- 1 2, OXYGEN AS PRESSURANT 
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TANK WALL TEMPERATURES, TANK CONFIGURATION D, 
TEST C003-10, HELIUM AS PRESSURANT 
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FIGURE 31b. ULLAGE PRESSURE AND PRESSURANT INLET TEMPERATURE 
HISTORIES, S-IV STAGE LOX TANK, HELIUM AS PRESSURANT 
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187260, NITROGEN AS PRESSURANT 
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FIGURE 35. COMPARISON OF PRESSURANT FLOWRATE PREDICTIONS BY 
TWO COMPUTER PROGRAMS WITH- EXPERIMENTAL RESULTS 
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BY TWO COMPUTER PROGRAMS WITH EXPERIMENTAL RESULTS 
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FIGURE 38. SCHEMATIC OF HEAT AND MASS TRANSFER 
CONDITIONS IN A PROPELLANT TANK 
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FIGURE 40. THE EFFECTS OF VARIOUS DESIGN PARAMETERS ON THE 
MEAN TEMPERATURE AT CUTOFF 
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